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1. Introduction 

The problem of the quantization of a many-body system attracted the attention of 
the physics community immediately in the first years after the introduction of the 
quantum mechanics [l]. Until today, different methods have been introduced to study 
this problem. Among the many we recall the Hartree-Fock method, the Thomas-Fermi 
approach, the Bohm-Madelung quantization, the second quantization by Heisenberg, 
the stochastic quantization by Nelson and others [2]. 

Concerning a non relativistic system, nonlinear Schrodinger equations (NSEs) are 
generally employed to take into account, through the nonlinearity, the interactions 
among the many particles. As a paradigm, the NSE with the cubic nonlinearity [U [5] 
has been widely used in literature to study, for instance, the Bose-Einstein condensation 
of alcali atoms hke ^Li, ^^Na and ®''Rb [6]. 

In [7], Bialynicki-Birula and Mycielski (BBM) introduced the nonlinearity —b ln(|^p) ip 
which was selected from the assumption of a separability condition among composed 
systems. Successively, in order to save only partially the superposition principle, valid 
in the linear theory, Weinberg [8] suggested a very general class of NSEs with an 
homogeneous nonlinearity. 

Motivated by these works, some experimental tests were carried out based on the 
neutrons interferometry [9l [ID] and ^Be"*" transition frequency spectroscopy [llj. All of 
them fail in pointing out the existence of nonlinear effects, suggesting in this way, that 
there is no real ground, at fundamental level, for such nonlinearities in the Schrodinger 
equation. 

Notwithstanding, NSEs can be usefully applied in the description of extended objects. 
In fact, depending on the form of the nonlinearity, some NSEs admit soliton solutions, 
i.e. non dispersive wave packets identifiable with quasi-classical particles. 
For instance, in the BBM equation a soliton solution named gausson has been studied 
[7[ [12]. This solution ([13] and references therein) represents the experimentally 
measurable charge or mass density p = l^/^p of an atomic nucleus, when it is considered 
as an elementary object (neglecting, in this way, its internal quarks structure). 
Remarkably, the BBM equation is closely connected with the Boltzmann-Gibbs (BG) 
entropy. In p^, it has been shown that the binding energy released in the splitting 
of an arbitrary wave function in non overlapping parts of the same form ip{x) 
"^■py^ ip(x—x.), with p. normalized according to p. = 1, is given by A E' = 6S'^*^(p), 
where S'^*^(p) = —J2iPi ^^Pt- This suggests that the BBM equation can be used 
in the thermodynamical description of spatially extended or self-interacting quantum 
mechanical objects. 

Also the BG entropy finds several applications in the study of quantum systems. For 
instance, in [HI [151 [IB] this entropy has been applied to generalize the Heisenberg 
uncertainty relation for a pair of conjugate observables. In [17] a worthwhile discussion 
concerning the BG entropy and the Fisher information measure has been advanced. 
Author clarifies the role of these entropies in the Smoluchowski diffusion process which 
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describes the underlying quantum kinetics of the hnear Schrodinger equation, 
Moreover, in recent years, an intensive research activity has been focused in the study 
of complex quantum systems which are intimately related to generalized versions of 
the BG entropy [HI [191 120], although the link between the microscopic properties of 
these systems and their statistical mechanics features is still an open question. In this 
respect, in [21] a microscopic derivation of a statistical system described by a general 
entropy has been considered. A microcanonical approach for a quantum system, at 
equilibrium, is studied in [221 ES] (and references therein) whereas, more recently in 
[21], author deals with the canonical quantization of a classical system governed by an 
arbitrary entropy whose kinetics is described by a very general nonlinear Fokker-Planck 
equation (NFPE). He obtained a wide family of NSEs where the structure of the complex 
nonlinearity is related to the form of the entropy of the ancestor classical system. The 
method has the drawback that, in consequence of the canonical quantization, this family 
of NSEs describes the time evolution of isolated quantum systems conserving energy and 
momentum. 

The purpose of the present paper is to accomplish the quantization of a many-body 
system of interacting particles in thermal contact with a reservoir (the background). 
This can be performed by using a generalized version of the stochastic quantization 
method initially proposed by Nelson [251 ISSl [27]. In this way, we derive a family of 
NSEs with a complex nonlinearity, describing damped and dissipative processes. We 
remark that complex nonlinearities arise whenever the quantum system obeys to a 
diffusive continuity equation. This, for example, occurs for the family of Schrodinger 
equations derived by Doebner and Goldin [28] as the most general class compatible with 
the linear Fokker-Planck equation, or, in [29j, where a NSE has been derived from a 
generalized Pauli exclusion-inclusion principle accounted through a suitable modification 
(in a nonlinear manner) of the particle quantum current. Also in [30], a class of NSEs 
with a complex nonlinearity has been obtained in the stochastic quantization framework, 
starting from the most general kinetics containing a nonlinear drift term and compatible 
with a linear diffusion term. Quite generally, all these NSEs can be transformed in other 
ones, with a purely real nonlinearity, by means of a gauge transformation of third kind 
[281 [331 [34]. 

Although sometime criticized [31], the Nelson method is, at present, argument 
of an intense research activity [35l [36l [371 [2H|- It merges together both the 
Newtonian mechanics and the diffusion process with the prospect to understand 
the quantum mechanics in a pure classical language. Typically, one assumes a 
quantum system undergoing an irreducible Brownian motion in the configuration space. 
Notwithstanding, the stochastic mechanics is not merely a mathematical description of 
a diffusion process. Two main aspects make the difference with the usual treatment of a 
standard diffusive process: the prescription associated to the reverse of the time axis in 
the diffusion process and the introduction of a stochastic acceleration, whose definition 
has been generalized in several papers [391 

In the following, we generalize the original Nelson method according to the following 
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two steps: 



(i) By assuming a nonlinear kinetic underling the stochastic diffusion motion, so that one 
can replace the forward and backward linear Fokker-Planck equations with a couple of 
nonlinear kinetic equations, in a manner compatible with the original NFPE describing 
the kinetic of a classical system; 

(ii) By generalizing the definition of the mean acceleration in the most general quadratic- 
form, compatible with a local expression of an Euler equation. 

In the next section 2, we recall the relationship among an arbitrary entropy and 
the corresponding NFPE, describing the kinetic of a classical system. In section 3, 
we start from this NFPE and derive the family of NSEs describing a quantum many- 
body system. In section 4, we investigate the time evolution of the main observables 
(Ehrenfest relations) whilst, in section 5, we inquire the existence of a special class 
of solutions, the solitons, which describe an isoentropic and non dissipative process. 
Finally, in the section 6, as an example, we discuss the quantization of a classical system 
described by the BG entropy. The conclusions are reported in the final section 7. 

2. Nonlinear kinetic equation 

In literature, there have been presented many different approach to relate, in a classical 
system, the production of entropy to a NFPE [32]. 

In the following, we consider a many-body classical system described with the field 
p{x,t), labeled by the vector x = [xi, X2, x^), and normalized in 



Quite generally, the time evolution of this field can by obtained, in the configuration 
space, starting from the continuity equation 




(2.1) 




(2.2) 



and assuming for the current the expression 
J = -pF{p) . 

The thermodynamic force F{p), defined through the relation 



(2.3) 




(2.4) 



is related to the functional 



U-DS{p) , 




where 




(2.6) 
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is the mean energy of the system, S{x) is the total energy for particle and D is the 
constant diffusion coefficient. 

Let us assume for the entropy S{p) the very general trace-form expression 

S{p) = - j dx j dp \nn{p) , (2.7) 

(throughout this paper we use units with the Boltzmann constant k-^ = 1). 
In equation (12.71) . the quantity K,{p) is a smooth function of the density field p{x, t). 
Depending on the expression of n{p), equation (12.71) encompasses different entropies 
already introduced in literature [H]. 

From (I2.2p - (I2.7I) we obtain the following continuity equation 

1^ + V ■ [p^drift -Dh{p) Vp] = , (2.8) 
where we have introduced the drift velocity 

Adrift = -V^(x) , (2.9) 
whilst the function /o(p) is given by 

/o(p)=p^ln/s:(p) . (2.10) 

Equation (12. 8p is a NFPE in the Smoluchowski picture since it describes the kinetics 
process in the configuration space rather than in the momentum space. 
We observe that in equation (12. Sp the total current J = Jdrift + «^difT is the sum of 
the standard linear drift current Jdrift = p*^ drift and the nonlinear diffusion current 
Jdiff = —Dfo{p) Vp which reduces to the standard Pick form Jpick = —D'Vp when 
k{p) = p. 

Since D is constant, equation (12.81) describes an isothermal diffusive process where the 
system is in thermal contact with a reservoir. In fact, it is known that for an isothermal 
process the diffusion coefficient plays the role of the inverse of a kind of temperature 
[32] . In this case, we easily recognize that equation (12. 5p defines a free energy for the 
system we are looking at. Actually , we can easily verify that the quantity C{p) is a 
Lyapunov function for the given problem since, according to equation (12.81) . we have 



dC 



V£{x)-DV 



6S 
~p 



2 



da; < , (2.11) 



where equality holds at equilibrium. 

It is worthy to recall that recently, different NFPEs have been related with 
generalized entropies of the kind (12. 7p by using several methods substantially equivalent 
to the one described above. In [12], generalized Kramers and Smoluchowski equations 
were derived by means of a variational principle, which minimizes the dissipative ratio of 
a generalized free energy equivalent to equation (12. 5p . The derivation of a Smoluchowski 
equation starting from a Kramers equation has been rigorously derived in [43] by mean 
of the Chapman-Enskog expansion in the high friction limit. Finally, in [H], a kinetic 
interaction principle has been proposed to derive a family of NFPEs starting from a 
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classical Markovian process describing the kinetics of a system of A^-body particles. The 
kinetic interaction principle is realized through a suitable factorization of the transition 
probability nit, x ^ y) = r(t, x, y) a{p) b{p') c{p, p') in terms of the population of the 
initial site x and the final site y, with p = p{t, x) and p' = p(t, y), where the functions 
a{p), b{p') and c{p, p') fix the nonlinearities in the NFPE and imposes the form of the 
entropy associated to the system. 



so that 



3. Stochastic quantization of the many-body system 

In the stochastic quantization framework, particles are assumed undergoing a time- 
asymmetric Markovian random process. The space coordinate x satisfies a couple of 
stochastic differential equations 

dx{t) = v^^\x,t)dt + dW^^\t) , (3.1) 

describing a forward (+) and backward (— ) process, respectively. They can be seen 
as a system of Langevin equations in the configuration space, where v^'^^x, t) are the 
forward and backward velocities, two vector valued functions depending on space and 
time, whilst the functions W^^\t) describe a Wiener process. 

According for equation (13. ip . x{t) is not a differentiable function. This means that 
there is no velocity as standard time derivative. Thus, one introduces a forward and a 
backward time derivative given by 

A) .^ + „(...V±Z,A. (3.2) 

v(^\x,t) = i^-j x{t). (3.3) 

Let us identify the field p{x, t) with the probability density of a quantum system at the 
position x{t). We assume that it fulfils the following forward and backward NFPEs 

|^ + V-[«WpT/^/(^Hp)Vp] = 0, (3.4) 

where, without lost of generality, we pose 

/W(p) = 2p^ln«:W(p) , (3.5) 

with K^^\p) two arbitrary functions. Clearly, the form of equations (13.41) is reminiscent 
of the expression (12.81) which describes the kinetic of a classical system. 

From equations (13. 4p . we see that forward and backward velocities are related each 
other in 

= -y(-) + V^(p) , (3.6) 

where 

g{p)=\n{n^^\p)K^-\p)) . (3.7) 
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V 



(+) 



V 



(-) 



(3.8) 



As customary, we introduce the mean velocity and the osmotic velocity according to the 
relations 

V = - fv'-^^ + v^~^^ , u = - 
2V / ' 2 

The former quantity represents the global velocity of the density shape of the quantum 

fluid, whilst the later, which takes the expression 

u = DVg{p), (3.9) 

has an intrinsic stochastic origin and is related to the spatial variation of the density. 
Endowed with the definitions fl3.8p . by averaging equations (13.41) . we obtain 



9p „ 



pv-Dfo{p) Vp 







a nonlinear continuity equation for the field p{x,t), where 



f-\p) 



p-^^'^^^ip) 



and 



/€(+)(p) 



(3.10) 



(3.11) 



(3.12) 



«:(-)(p) • 

Equation (13.101) is formally equivalent to the NFPE (12.81) . According to this 
correspondence, the expression of the function k{p) is fixed through the entropy of 
the ancestor classical system although the same entropy does not determine univocally 
the functions K^^\p). From a mathematical point of view, this means that there are 
infinitely possible choices for the functions K^'^^p), all compatible with the form of the 
entropy (12. 7p . according to equation (13.121) . Note that equation (I3.10p is consistent with 
the kinetic interaction principle proposed in [H], where a very general NFPE, whose 
equation (I3.10p is a particular case with the linear drift, has been derived starting from 
a master Pauli equation describing a Markovian process. 

Let us now introduce the stochastic acceleration. In the original Nelson work it is 
defined in 



a 



1 
2 



(+) 



(-) 



V 



(+) 



(3.13) 



i.e. by means of the algebraic mean of the quadratic operators (d / dt)^^\d / dt)^ ^ and 
{d/dtY-\d/dtY+^ acting on the vector x{t). The non uniqueness of the definition of 
a has been noted in [39] and alternative definitions have been proposed in [391 HO]- 
Hereinafter, we assume the most general quadratic form, in the forward and backward 
time derivatives, compatible with a local expression for an Euler equation describing 
the dynamics of v. Thus, we pose 



a 



(+) 



-) 

dt J 



(+) 



V 



(-) 



(-) 



V 



{+) 



+ V 



(-) 



V 



(3.14) 



with the condition t + 2 p + v = 1. 

We observe that the original Nelson definition is recovered from equation (I3.14p by 



Stochastic quantization of interacting classical particles system 



8 



posing T = u = and /i = 1/2 whilst, the definition advanced in is obtained by 
posing T = u = 1/8 and /i = 3/8. Moreover, in the framework of a stochastic Lagrangian 
formulation [40], compatible with the variational stochastic calculus introduced in |46j . 
another definition of a has been advanced. It follows from equation fl3.14p by posing 
T = u = 1/2 and fi = 0. 

As customary, let us introduce the potential field S(a3,t) related to the mean 
velocity through the relation 



m V 



(3.15) 



Remarkably, this equation is formally equivalent to the relation (12. 9p among the classical 
drift velocity and the particle energy spectrum although the fields S(a3,t) and S{x) 
represent two different physical quantities. In this sense, the quantum mean velocity v 
assumes a different meaning, in the quantum picture, with respect to the meaning that 
the classical drift velocity Vdrift has in the classical picture. 

Taking into account that the both fields v and u are given through the gradient of the 
scalar functions g{p) and S(a;,t), respectively, and recalling the Newton's law, from the 
definition fl3.14p one obtain 



dv „ du „, 

— + V -Vv + ai — + ai V{u 
at at 



V 



+ DaiAv 



+ ■ Vii + D a2 Au 



1 



VV = 



m 



(3.16) 



where ai = t — u and a2 = t — 2 fi + u. 

In effects, (13.161) is an Euler equation for the mean velocity field v{x,t). We remark 
that, since the osmotic velocity is a function of the field p{x, t), the quantity du/dt can 
be easily handled by using the continuity equation (I3.10p . 

By employing the definitions (13.91) and (I3.15p . the iiJu/er equation takes the expression 



'dt 



+ 



where 



h{p) 

f2{p) 



(VS) 

2m 



aiD 1 - p 



+ /i(p) AS + Mp) {Vpf + /3(p) A p + y = 



d9{p) \ 
dp ) 



(3.17) 



ai 



dfojp) dgjp) a2 f dg{p) \\ d'^gjp) 
— z z \ I — — I + a2 - . 



h{p)=mD^ 



dp 
fli /o(p) 



0-2 



dp ' 2 

dgjp) 
dp 



\ dp J 



dp^ 



(3.18) 



As known, in the hydrodynamic representation of a quantum system [HI HS], the 
quantum effects are take into account through the introduction of the quantum potential 
[/q(p) = — (/l^/2 m) A^yp/^yp. In the original Nelson formulation, the potential t/q(p) 
arises from the particular expression assumed by the osmotic velocity that, in that case, 
is given by 

u = D\Yip. (3.19) 
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In the nonlinear case studied in this paper, the same form fl3.19p can be recovered by 
posing 



kW(p) = y^Mp) , '^^-\P) = , (3.20) 

although, according with the non uniqueness of K^^\p), the choice (KW\f is not unique. 
In this way, equation fl3.17p simplifies in 

^ + + U,{p) + W{p) + V = 0, (3.21) 

where 

W{p) =mD'Jv- (/o(p) Vp) , (3.22) 
and 02 = —{h/2mDy. 

Equations f l3.10p and f l3.2ip govern the time evolution of the fields p{x,t) and T,{x,t). 
Actually, fl3.10p is a quantum continuity equation for the field p{x, t) describing the 
kinetics of the quantum fluid a la Bohm-Madelung, where the quantum current, 
according to the position fl3.15p . assumes the expression 

j = p^-Dfoip)Vp. (3.23) 
m 

Differently, fl3.2ip can be interpreted as a nonlinear Hamilton- J acobi equation for the 
potential field S(a;, t) whose dynamics is ruled through the quantum potential, the extra 
nonlinear potential W{p) and the external potential V . 
Finally, with the position 

iP{x, t) = ^/p{x, t) exp (^^ E{x, t)^ , (3.24) 

both equations fl3.10p and f l3.2ip can be joined to obtain the following NSE 

th^ = R^, (3.25) 
where H is the non-Hermitian operator 

R = -^A + Aip) + V , (3.26) 
2 m 

with a nonlinearity A(p) depending only on the field p{x,t). It is given by 

A(p) = ^V-(/o(p)Vp) , (3.27) 

where the complex coupling constant takes the expression 

h 

X = mD^ai+i-D. (3.28) 

The family of NSEs fl3.25p describes, in the mean filed approximation, the dynamics 
of a quantum system of identical particles interacting with a thermal reservoir. The 
nonlinearity A(p) is fixed through the entropy ( 12. 7p which governs the kinetics of the 
ancestor classical system and consequently determines also the continuity equation of 
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the quantum system. 

In conclusion, let us make a few considerations: 

- First, the original Nelson derivation is recovered by posing t = u = and fi = 1/2. 
This imply 02 = —1 and we recover the well know relation D = h/2m among the 
diffusion coefficient and the mass of the particles of the system [25]. Moreover, according 
to the original proposal, we must pose K^+^(p) = K^^\p) = ^p. As a consequence /o(p) 
vanishes and equation f l3.25p reduces to the ordinary linear Schrodinger equation. 

- Second, when ai vanishes, equation fl3.25p describes a family of NSEs with a 
purely imaginary nonlinearity. In this case only the continuity equation acquires an 
extra nonlinear term (the diffusive one) whilst the Euler equation maintains the same 
expression of the Bohm-Madelung theory. 

- Finally, by means of a nonlinear gauge transformation [28l [331 El! > family of 
NSEs (13.251) can be transformed in another one with a purely real nonlinearity. 

In fact, by performing the transformation 

/o(p) 



P 



dp 



ip{x, t) — > (f){x, t) = ip{x, t) exp (y^j^ J 

the nonlinearity f l3.27p changes according to 

A(p) ^ A(p, a) = (p) Aa + g^ip) (Vpf + g^{p) Ap 



(3.29) 



(3.30) 



with 



g,{p) = -D foip) 



g^[p) =mD' 



1 ( foip) 

2 V P 



+ 



^ dfojp) 
P dp 



(3.31) 



^3(p) =mD^ai 



Up) 



In equation fl3.30l) . a{x,t) is the phase of the new field (j){x,t), related to the phase of 
the old field il){x^t) through the relation 



a{x,t) = T.ix,t)-mD [ dp . 

J P 



(3.32) 



We observe that the new nonlinearity A(p, a), albeit purely real, now depends through 
the both fields p{x,t) and cr{x,t). 

4. Time evolution of observables 



The family of NSEs fl3.25p describes dissipative and damped processes. To show this we 
derive the Ehrenfest relations for the main observables of the system. 
In passing, we notice that, as a consequence of the continuity equation fl3.10p . the 
normalization of the wave function 



p{x, t) dx 



(4.1) 
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is preserved in time. We recall that nonlinear evolution equations, in general, are not 
ray-invariants. This means that, given a solution of equation (13.251) . Aip with A a 
constant, does not fulfill the same equation. Therefore, in this case the normalization 
condition (14.11) can be accomplished by fixing appropriately one of the free parameters 
of the model. As discussed in ||8j, the ray-invariance can be restored when equation 
( I3.25P is homogeneous in p{x, t). 

Let us now introduce the linear momentum P = j p VS dx. By using equations 
( I3.10p and (I3.2ip . we can derive the following continuity equation 

^ + V ■ T = 7?-p , (4.2) 
for the density p = p VS, where the momentum-stress tensor T has entries 



T 



L d.T. d,T. + d,p d,p + G'i(p) 

m Imp 



with 



Gi{p) 



4 m 



Ap + aimD'^ p — 



fpT- -V-(/o(p)Vp 

J dp Ip \ 



dp . 



(4.3) 



(4.4) 



Clearly, the momentum source term IZp is responsible of the non conservation of P. Its 
expression is given by 

1Zf = G2{p)VJ:- pVV , (4.5) 

with 



G2(p) = Z^V-(/o(p)Vp) . 



(4.6) 



In the same way, by defining the energy of the system E = — J p (d'E/dt) dx, we 
can derive the evolution equation for the density e = —pdT,/dt, which reads 
de 



dt 



p7^I 



where the energy current Je is given by 



dt 



m ot Am p dt 
The energy source term takes the form 



^ ^ , , 9S ^ , .dp dV 



with 



d 

Gsip) = aimD'^ p — 
dp 



iv. 

.p 



Up) 
p 



Vp 



(4.7) 



(4.8) 



(4.9) 



(4.10) 



and the quantities dp/dt and dT,/dt are computed through equations (13.101) and (I3.2ip . 
respectively. 

By integrating equations (14. 2 p and (14.70 on the whole configuration space and 
assuming uniform boundary conditions on the field p{x,t) we obtain the relations 
dP I \ dE 

-dT - ' -d^ 



(4.11) 
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(4.12) 



where (T^x) = / TZ^dx. 

A further couple of relations easily derivable from equations (14.21) and (13.101) is given by 

dL / 7^ \ dx^c P 

dt \ ^ / dt m ' 

with L = j X X pdx the total angular momentum and x^c = J pxdx the mass center 
of the system. 

Thus, as stated by relations (I4.11I) - (I4.12I) . equation (13.251) describes a quantum 
dissipative process due to the presence of the sources TZp and TZe- Notwithstanding, 
depending on the form of /o(p), there could exist special classes of solutions describing 
conservative systems. 

5. Soliton solutions 



Among the many solutions of a NSE, solitons are of particular interest in Physics. 
As known, soliton solutions arise when the dispersive effects, due to the term Aip are 
balanced by the nonlinear term A(p). In this case, the many particle system moves 
coherently giving origin to an extended object with a particle-like behavior preserving 
its shape in time. 

When the system is constrained in one dimensional space, without the external 
potential V = 0, a soliton solution can be obtained through the following ansatz 



exp 



mv'^t 



(5.1) 



_hV'^' 2 
where ^ = x — vt and v is the soliton velocity. 

Accounting for the relations d/dt = —vd/d^ and d/dx = d/d^, equations (13.101) and 
(I3.2ip form a couple of ordinary differential equations for the scalar fields s(^) and p(^) 



ds 



dp 



P ~T7 ~ ''^ ^ fo{p) -77 ~ P = ^ ) 



ds 
d^ 



d^ 



H — mv 
2 



By introducing the function 

\P 

we can obtain the following first order differential equation 

dy{p) 



dp 



where 



Flip) = - 
P 



F,{p)y{p) + F2{p) = , 

/o(p) (/o(p) + 2ai) - ao + 2 ai p (dfoip) / dp) 



ai /o(p) - ao 



F2{p) 



1 2 iv/D) 



pa-i fo{p) - ciQ ' 



(5.2) 
(5.3) 

(5.4) 
(5.5) 

(5.6) 
(5.7) 
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and Oo = {h/2m D)'^. 

Equation (15.51) is solved by quadrature so that, accounting for equation (15. 4p . we get 

-1/2 

dp , (5.^ 



^ = 1 i Ce-/^^^"')'^''' - I F^ip') e/^^(''")'^^" dp' 



P L 

where C is the integration constant. 

This equation defines imphcitly the shape of the sohton. Physically relevant solutions 
arise when p(^) vanishes in a sufficiently fast way for C, —>■ ±cxd, so that the integral 
I-^ PiO d-^ converges. Clearly, this last condition selects the form of /o(p) and 
consequently the possible NSEs admitting soliton solutions. 

We remark that this special class of solutions is not dissipative nor damped. They are 
steady states for the system under inspection. In fact, we easily verify that 

+ 00 +00 

— oo — oo 

and according to the boundary condition p(±oo, t) ^ we obtain that 
In the same way 

+ 00 +00 

/" (is , f ds ,^ 1 2 1 

h = / p —— dx = —V / p — dt mv , (5.11) 

J dt J ^ d^ ^ 2 ' ^ ' 

— oo — oo 

so that 

Actually, the solution (15. 8p describes a particle-like object traveling with constant 
velocity v, momentum P = mv and energy E = mv'^/2. It is also easily to verify 
that this solution describes an isoentropic process. This is consistent with the meaning 
of the entropy as a measure of the information. In fact, since the soliton does not change 
its shape in time, it is reasonable that the information (entropy) contained in it does 
not change in time, accordingly within the relation dS{p)/dt = 0. 

6. An example 

As an example, let us consider the quantization of a classical system described by the 
Boltzmann-Gibbs entropy 

S^^{p) = - j plnpdx , (6.1) 

which is obtainable from equation ( 12. 7p by posing k{p) = p/e. This implies /o(p) = 1, 
so that the kinetics of the classical system is governed by the linear Fokker-Planck 
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equation. 

The quantum system is described by the following homogeneous NSE 

^h^ = -^Ai, + \^i; + V^P, (6.2) 
ot 2 m p 

where A is given in equation ( 13.28^ . This evolution equation belongs to the family 
of NSEs study in [28] and derived on the physical ground as the most general class 
of nonlinear quantum evolution equations compatible with the linear Fokker-Planck 
equation. 

By employing the gauge transformation (I3.29p . equation (16.21) reduces to a NSE with a 
purely real nonlinearity 

ik^ = -—A(t) + W{p,(j)(t) + V(t), (6.3) 



where 



and 



dt 2 m 



W{p,a) = -DAa + -mD^ [ ] +aimD^^, (6.4) 



1 n2^VpV , . ^n2Ap 



(7 = 11 — mD In p. (6.5) 

By posing ai = (ao — l)/2 and V = 0, equation (16. 2p admits a physical solution whose 
density and phase read 



p{x,t) = ^exp(^-^mc, (I) ^ 



(6.6) 



S(x,t) =C2 + J 



^""^ {Am^D'' + h^) + -mx' 



Im^^d)', (6.7) 



16 ' '2 

respectively, where A, ci > and C2 are constants. Accounting for the homogeneity 
of equation (16. 2p we can accomplish the normalization of p{x, t) by posing A = 
(mci/4 7r)^/^. 

Solution (I6.6p - (I6.7I) describes a dispersive wave function spreading out for t — > +oo, 
since p{x, +oo) — >• 0. Correspondingly, we have also P{+oo) — > and E{+oo) 0, 
indicating that this solution describes a not reversible process, damped and dissipative. 
It is also possible to verify that the corresponding entropy is given by 

S^^ = \nt + Ct-lnA , (6.8) 

with C = 12 (7r/m Ci)^/^, which grows in time. 

With the same choice for the constant ai and for the potential V, equation ( 16.21) 
admits a soliton solution. Its shape has the form 

p(e) = Ae-f«' , (6.9) 



where a = {AmvY / {Am^ D"^ + h^) and A = yJajA'K. The phase of the wave function 
is given by equation (15. 2p and takes the expression 

1 1 
Y.{i,t) = —^mDai^ ^mvi- -mv^t . (6.10) 
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The Gaussian shape (16.91) is known in hterature as a gausson [TJ [T2] and it was derived 
firstly as the sohton solution of the BBM equation, although, in that case, the phase 
assumes a different expression. Solution fl6.9l) -( l6TT0l) describes a not dissipative 
reversible process with a constant entropy 



7. Conclusions 

In the stochastic quantization framework, we have approached the problem of the 
quantization of a many-body classical system whose kinetics is described by a general 
NFPE compatible with an arbitrary trace form entropy. By a suitable modification of the 
original Nelson method we have derived a family of NSEs describing, in the mean field 
approximation, a system of collectively interacting particles. The time evolution of the 
mean field ip{x, t) is given by means of a non-Hermitian operator H whose nonlinearity, 
describing dissipative and damped effects, takes into account the collective interactions 
among the many particles and the background. This is a substantial different situation 
respect to the canonical approach studied in [21] which deals only with conservative 
systems. 

We have investigated the existence of a special class of solutions representing a localized 
object with particle-like properties originated by the collective motion of the many 
particles of the system. This solution describes a reversible and non dissipative process 
without production of entropy. 

As an explicit example, we have studied the quantization of a classical system governed 
by the BG entropy obtaining, in this way, a NSE with a complex nonlinearity. This 
example emphasizes the difference among the quantization method employed in the 
present paper and the canonical one used in [21]. There, starting from a classical system 
governed by the BG entropy, the resulting quantum evolution equation is given by the 
linear Schrodinger equation which describes the dynamics of an isolated system. This is 
consistent with the widely accepted opinion that the BG entropy governs systems with 
weakly and local interactions. Differently, in the present case, although the particles in 
the system must still be considered weakly interacting (due to the presence of the BG 
entropy), there is a further interaction between the system and the background which 
is accounted for the nonlinearity in the evolution equation. 
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